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A pedagogically instructive experimental procedure is sug- 
gested for distinguishing between different damping terms in a 
weakly damped oscillator, which highlights the connection be- 
tween non-linear damping and initial-amplitude dependence. 
The most common damping terms such as contact friction, air 
resistance, viscous drag, and electromagnetic damping have 
velocity dependences of the form constant, v, or v 2 . The cor- 
responding energy dependences of the form y/E, E, or E\f~E 
in the energy loss equation give rise to characteristic depen- 
dence of the amplitude decay slope on the initial amplitude. 



I. INTRODUCTION 

The most commonly studied source of damping in os- 
cillating systems is a resistive force F ros = — mv/r pro- 
portional to velocity, such as due to fluid drag at very 
low relative speed, when the Reynolds number R is of 
the order 1 or less. The equation of motion involves a 
linear differential equation of the type 



H h UJnX = 



(1) 



which yields an exponential decay xo(t)/xo{0) — e~*/ 2r 
for the normalized amplitude of oscillations. The decay 
is independent of the initial amplitude, a characteristic 
of the linear nature of the system. Here x(i) refers to 
the instantaneous displacement of a block in a spring- 
mass oscillator or the angular displacement of a physical 
pendulum. 

For weakly damped oscillations, it is convenient to con- 
sider the oscillator energy E averaged over one cycle, 
during which period the amplitude xo{t) is nearly con- 
stant. With x(t) — xosmujt, we have E — 17110^x^/2, 
and the change in the oscillator energy over one cycle 
is AE — J F lcs .vdt = -E T/t. This yields the linear 
energy loss equation 



dE 
~dt 



E 

T 



(2) 



over long time scales compared to the oscillation time 
period. The normalized energy E{t)/E{Q) — e~ l l T also 
decays exponentially, independently of the initial energy. 

In most laboratory oscillators, however, non-linear 
terms actually dominate in the energy loss equation. A 
\f~E~ term results from the nearly velocity-independent 
dry friction, for instance at the pivot of a physical pen- 
dulum or due to a brake pad, whereas air resistance leads 
to a E\[E term as the Reynolds number is typically of 
the order of 1000, in which regime the air resistance is 



proportional to v 2 . Another potential contribution to the 
energy loss equation arises from a centripetal correction 
to the normal reaction at the pivot, and therefore to the 
friction. As the average centripetal force F c cx w^ax x E, 
this leads to a E\f~E correction to the contact friction 
contribution in the energy loss equation. 

In certain situations, even electromagnetic damping in 
an oscillating system leads to a \J~E damping term, as re- 
cently reported for a simple experimental setup involving 
a magnet oscillating through a coil. [Q Similar damping 
behavior is expected in an oscillating pendulum consist- 
ing of a copper or aluminum disk, periodically passing be- 
tween the pole pieces of two strong magnets placed near 
the mean position, provided the conducting disk passes 
through the magnetic field in a very short time compared 
to the oscillation time period. On the other hand, a con- 
tinuously acting electromagnetic damping force, as on 
a conducting spinning wheel placed between two strong 
magnets, leads to a E damping term. 

In practice, the oscillator energy loss is therefore typi- 
cally described by the following non-linear equation 



^E = -k^E-k'E-k"EVE 
at 



(3) 



In this article we describe a quantitative study of the ex- 
pected initial-condition dependence of the amplitude de- 
cay due to the non-linear damping, and discuss a practi- 
cal application of this initial-condition dependence. Even 
though the different damping terms have characteristic 
decay signatures — a linear decay of amplitude with time 
for \f~E~ damping, an exponential decay for E damping, 
and an inverse power decay for E\f~E damping — a prac- 
tical difficulty often encountered is how to distinguish 
between several weak damping terms present simultane- 
ously. Towards this end, we suggest a sensitive experi- 
mental procedure for quantitatively identifying the sepa- 
rate damping contributions in Eq. (3). Furthermore, as 
this procedure requires only the initial amplitude decay 
and not the full decay over long time, it is especially use- 
ful in situations where very small oscillation amplitudes 
cannot be reliably obtained. 

Either separately, or in combination, the different 
damping terms have been considered in several earlier 
studies. These include a harmonic oscillator with sliding 
friction and viscous force, || an oscillating sphere 
with fluid drag, [pi] and a physical pendulum with air 



resistance 
damping. |11 
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FIG. 1. A magnet (NS) attached to an oscillating system 
passes through a coil (C) periodically, generating a series of 
emf pulses. 



II. EXPERIMENTAL SETUP 

We have used a common laboratory setup for studying 
electromagnetic induction to monitor the oscillation am- 
plitude. An oscillating magnet passes through a coil peri- 
odically, generating a series of electromotive force (emf) 
pulses. In a practical realization of this concept, a rigid 
semi-circular frame of aluminum, pivoted at the center 
(O) of the circle (see Fig. 1) and free to oscillate in 
its own plane about a horizontal axis through O, has a 
rectangular bar magnet mounted at the center of the arc 
passing through a coil C of suitable area of cross section. 
A convenient way of monitoring the induced emf pulses, 
and therefore the oscillation amplitude, is through a PC 
interface, which can be readily realized by low-cost, con- 
venient data-acquisition modules available in the market. 
The amplitude can also be directly monitored by con- 
necting the oscillator shaft to a precision potentiometer, 
and recording the instantaneous voltages through a PC 
interface. |l(J 

The un- 
derlying electromagnetic induction phenomenon in this 
oscillating system has been discussed earlier. [[[] The in- 
duced emf is significant only in a very narrow angular 
range about the mean position, and if the angular am- 
plitude is not too small (9 ^> 5°), the angular velocity 
of the bar magnet is very nearly constant in this narrow 
angular range. The peak emf fmax is then approximately 
given by 



~d9 



2tt 
— f 
T 



(4) 



where T is the time period of (small) oscillations, and the 
maximum angular velocity w max = (Att/T) sin(0o/2) « 
(2w/T)9o, if 9q/2 (in radians) is small compared to 1. 
Thus, the peak emf £ max provides a measure of the angu- 
lar amplitude 6q, & n d the oscillator energy E 
is directly proportional to 9q. 



1 T, ,2 

2 ^max 



III. DECAY OF OSCILLATOR ENERGY AND 
AMPLITUDE 

While a solution of Eq. (3) is easily obtained, it is 
more instructive to proceed in steps, and we consider 
three separate cases. 



A. k' = k" = 

When only a \f~E~ damping term is present, as due to 
contact friction and intermittent electromagnetic damp- 
ing, the solution of Eq. (3) yields a linear decay of the 
normalized oscillator amplitude with time 



%(t) _ y/E® 



9o(0) y/E$) 



= l-at 



(5) 



where the decay slope a = k/2y/E(0) is inversely pro- 
portional to the initial amplitude. A phase-space plot 
between the normalized angular momentum L/Lo and 
angular displacement 9/0q is a linear spiral, and the num- 
ber of cycles executed before the oscillations come to a 
stop is proportional to the initial amplitude. 



B. k" = 

When a resistive force proportional to velocity is in- 
cluded, which would be appropriate for low relative ve- 
locities when the Reynolds number R ~ 1, the solution 
is modified to 



k + k'^/E(0) 



which leads to an exponential decay E(t)/E(0) = e 
in the limit k — *■ 0. It is instructive to consider the k' - 
limit. Expanding the exponential term, we obtain 



%{t) _ y/Elfi 



9 (0) y/WO) 



l-at + bt 2 + 



(6) 

-fe'i 

-> 



(7) 



where a = k/2 y / E(0) + k'/2 and b/a = k'/A. The E 
damping term contributes to the deviation from linearity, 
and as k' — > 0, all second- and higher-order terms vanish, 
leaving only the linear term, as in Eq. (5). 



C. General case 

We again consider the small-time behaviour of the os- 
cillator energy E{t) or amplitude &o(t) given by Eq. (7), 
and substituting in Eq. (3) we obtain 
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FIG. 2. The normalized peak voltage V r max (t)/l / max(0) vs. 
time t for the short-circuit (+) and open-circuit ( x ) cases. 

2a=-^= + k' + k"^E(tij , (8) 

b k> fcVgo) 

a ~ J + 2 ' (9) 

As expected, the initial energy E(0) appears with both 
the non- linear terms k and k", and hence both a and b 
depend on the initial amplitude. More importantly, the 
decay slope a depends on the damping constants k, k', 
and k" , which implies that the linear part of the am- 
plitude decay itself contains information about all the 
damping terms. 

A typical amplitude decay with time is shown in Fig. 
2 for both the open- and short-circuit configurations. A 
dominantly linear decay is seen in both cases, with a 
deviation from linearity becoming pronounced at large 
times. In the open-circuit configuration, the damping 
is due to friction and air resistance. When the coil is 
short-circuited (through a low resistance (220 f2) so that 
the emf pulses can be monitored by tapping the voltage 
across the resistor), intermittent electromagnetic damp- 
ing is activated due to the induced eddy current when 
the magnet passes through the coil. 

Least-square fits with the quadratic form given in Eq. 
(7) yield a = 49 x 10" 4 sec -1 and b = 51 x 1CT 7 sec~ 2 for 
the open-circuit case, whereas for the short-circuit case 
a = 117 x 1(T 4 sec" 1 and b = 214 x 1CT 7 sec~ 2 . The 
linear decay rate a is significantly larger in the short- 
circuit case due to the additional contribution to k from 
electromagnetic damping. From Eq. (9), we note that 
the ratio b/a depends only on the air-resistance damp- 
ing coefficients k' and k" . If the electromagnetic damp- 
ing only modifies the k term, Q then this ratio should 
be indentical for both cases, provided the initial ampli- 
tude is identical. The near doubling of this ratio in the 
short-circuit case is therefore a clear indication that elec- 
tromagnetic damping also introduces small E and E\[E 
damping terms which modify the quadratic coefficient b. 



FIG. 3. The decay of the peak emf with time for different 
initial angular amplitudes shows no change in the decay slope, 
hence indicating no observable effect of air resistance. 

IV. INITIAL-AMPLITUDE DEPENDENCE OF 
THE DECAY SLOPE 

Equation (9) shows that both k' and k" contribute to 
the coefficient b of the quadratic term in Eq. (7). There- 
fore, it is not possible to distinguish between the E and 
E\^E damping terms from the initial amplitude decay, 
unless the dependence on the initial energy (amplitude) 
is taken into account. The form of Eq. (8) suggests 
a pedagogically instructive experimental procedure for 
identifying the signature of the weak damping terms k, 
k' , and k" by studying the initial amplitude dependence. 

If the amplitude decay is written as 9o(t) = #o(Q) — 
at + /3t 2 , and the peak emf £ max (i) = K9 (t) as 

£nax(<) = £max(0) - at + (3t 2 , (10) 

then from Eqs. (7) and (8), the slope a — K8 (0)a of 
the linear part of the decay can be written as 

a ^^[k + k' ce Q (o) + k" C 2 6 2 (0)] , (11) 



where y/E(0) = C6* o (0) defines the constant C. There- 
fore, a constant, linear, and quadratic dependence of the 
slope a on the initial amplitude 0q(0) will be signatures 
of the weak damping terms k, kl and k" , respectively. 

Figure 3 shows the amplitude decay for different initial 
amplitudes 9q(0) in the open-circuit configuration. There 
is no noticeable change in the slope, as confirmed from 
the plot of the decay slope a vs. 9o(0) shown in Fig. 4. 
The slope remains constant at 5.15 x 10~ 2 V/sec, and 
this rules out both the damping terms k' and k" . It thus 
appears that air resistance and centripetal correction to 
friction do not have any observable contribution to the 
amplitude decay, and the friction contribution is given 
by (K/2C)k = 5.15 x 10~ 2 V/sec. 
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FIG. 4. The decay slope a shows no change with the initial 
angular amplitude #o(0), showing no observable contribution 
to the amplitude decay from air resistance. 



FIG. 6. Increase in the decay slope a with the initial 
amplitude #o(0) is an indicator of the damping terms k' and 
k". 
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FIG. 5. A variation in the decay slope with the initial an- 
gular amplitudes is clearly observed when a 8| x 11" card- 
board sheet is attached to the oscillator to provide air resis- 
tance. 

When a 8^" x 11" cardboard sheet was attached to 
the oscillator frame to provide air resistance, the am- 
plitude decay showed a clear dependence on the initial 
amplitude 6>o(0), as seen in Fig. 5. Least-square fit 
of the peak-emf decay with the quadratic form of Eq. 

(10) yields the decay slope a, which shows a signifi- 
cant increase with 6*0(0) (see Fig. 6). Least-square fit 
of the decay slope a with the quadratic form of Eq. 

(11) yields (K/2C)k'C = 3.1 x 10~ 2 V/scc/deg and 
(K/2C)k"C 2 = 3.2 x 10- 3 V/sec/deg 2 for the linear 
and quadratic coefficients, respectively. It is instructive 
to note that the air resistance contribution to the decay 
slope a due to the cardboard sheet is comparable in mag- 
nitude to that due to friction at the pivot. The presence 
of a small k' term, although ideally this term should van- 
ish if the air resistance were purely proportional to v 2 , 
has also been noted earlier. |fl(l 



V. SUMMARY 

Qualitatively different damping terms can be distin- 
guished from each other by studying the variation of 
the amplitude decay slope with the initial amplitude. A 
constant, linear, and quadratic dependence of the decay 
slope on the initial amplitude are signatures of the E 1 / 2 , 
E, and E 3 ^ 2 damping terms, respectively. A quantitative 
determination of the contributions of contact friction and 
air resistance due to an attached vane is demonstrated. 
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